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In these Lectures I review possible constraints on particle physics models, ob- 
tained by means of combining the results of collider measurements with astro- 
physical data. I emphasize the theoretical-model dependence of these results. I 
discuss supcrsymmetric dark matter constraints at colliders (mainly LHC) in var- 
ious theoretical contexts: the standard Cosmological-Constant-Cold-Dark-Matter 
(ACDM) model, (super)string-inspired ones and non-equilibrium relaxation dark 
energy models. I then investigate the capability of LHC measurements in assert- 
ing whether supcrsymmetric matter (if discovered) constitutes part, or all, of the 
astrophysical dark matter. I also discuss prospects for improving the constraints 
in future precision facilities, such as the International Linear Collider. 



1. Introduction 

In the past decade we have witnessed spectacular progress in precision 
measurements in astrophysics as a result of significant improvements in 
terrestrial and extraterrestrial instrumentation. From the point of view of 
interest to particle physics, the most spectacular claims from astrophysics 
came in 1999 from the study of distant (redshifts z ~ 1) supernovae (SNe) 
of type la by two independent groups . These observations pointed to- 
wards a current era acceleration of our Universe, something that could 
be explained either by a non-zero cosmological constant in a Friedman- 
Robertson- Walker-Einstein Universe, or in general by a non-zero dark en- 
ergy component, which could even be relaxing to zero (the data are con- 
sistent with this possibility). In the past five years 2 many more distant 
(z > 1) supernovae have been discovered, exhibiting similar features as the 
previous measurements, thereby supporting the geometric interpretation 
of the acceleration of the Universe today, and arguing against the nuclear 
physics or intergalactic dust effects. 

Moreover, there is strong additional evidence from quite different in ori- 
gin and thus independent astrophysical observations, those associated with 
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Figure 1. The energy content of our Universe as obtained by fitting data of WMAP 
satellite. The chart is in perfect agreement with earlier claims made by direct mea- 
surements of a current era acceleration of the Universe from distant supernovae type la 
(courtesy of http://map.gsfc.nasa.gov/). 

the WMAP measurements of the cosmic microwave background radiation 
(CMB) 3 , as well as baryon acoustic oscillations measurements 4 . After 
three years of running, WMAP measured CMB anisotropies to an unprece- 
dented accuracy of billionth of a Kelvin degree, thereby correcting previ- 
ous measurements by the Cosmic Background Explorer (COBE) satellite 5 
by several orders of magnitude. This new satellite experiment, therefore, 
opened up a new era for astro-particle physics, given that such accuracies 
allow for a determination (using best fit models of the Universe) of cosmo- 
logical parameters 6 , and in particular cosmological densities, which, as we 
shall discuss in this review, is quite relevant for constraining models of par- 
ticle physics to a significant degree. All these measurements point towards 
the fact that (more than) 73 % of the Universe vacuum energy consists of 
a dark (unknown) energy substance, in agreement with the (preliminary) 
supernovae observations (see fig. 1). This claim, if true, could revolutionize 
our understanding of the basic physics governing fundamental interactions 
in Nature. Indeed, only a few years ago, particle theorists were trying to 
identify (alas in vain!) an exact symmetry of nature that could set the cos- 
mological constant (or more generally the vacuum energy) to zero. Now, as- 
trophysical observations point to the contrary. The WMAP satellite exper- 
iment determined the most important cosmological parameters that could 
be of relevance to particle physicists, namely 6 : the Hubble constant, and 
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thus the age of the Universe, the thickness of the last scattering surface, the 
dark energy and dark matter content of the Universe (to an unprecedented 
accuracy) (c.f. figure 1), confirming the earlier claims from supernovae la 
data and provided evidence for early re-ionization (z ~ 20), which, at 
least from the point of view of large scale structure formation, excludes 
Warm Dark Matter particle theory models. Its measurements have been 
independently confirmed recently by Baryon Oscillation measurements 4 . 

In this review I shall first describe briefly the above-mentioned astro- 
physical measurements, and then use them to constrain 7 some particle 
physics supersymmetric models, in particular (i) the minimally supersym- 
metric model, constrained by its embedding in a minimal supergravity 
model (mSUGRA) 8 , (ii) a class of string-inspired cosmologies 9 , with non- 
trivial dilaton and/or moduli fields, and (iii) a class of stringy relaxation 
dark energy models 10 . I shall give a critical discussion on the derived con- 
straints, and their strong theoretical-model dependence, and discuss the 
capability of observing Supersymmetry in colliders after the latest astro- 
physical data. In this respect I will also discuss the importance of the — 2 
experiments of the muon gyromagnctic ratio, in light of some very recent 
(2006) results 11 , pointing towards a discrepancy between the measured 

— 2 value from the one calculated within the standard model. 



2. Astrophysical Measurements and Facts 

On Large scales our Universe looks isotropic and homogeneous. A good 
formal description, which does not depend on the detailed underlying mi- 
croscopic model, is provided by the Robertson- Walker (RW) metric, ac- 
cording to which the geometry of the Universe is described by means of the 
following space-time invariant element 12 : 



ds 2 = -dt 2 + a(t) 2 Rl 



dr 
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(2.1) 



where a(t) = S^- = is the scale factor, H = | is the Hubble Pa- 
rameter, t is the Cosmological Observer time, Ro denotes the present-day 
scale factor, z — is the rcdshift, k denotes the Spatial Curvature, which (by 
normalization) can take on the values: k=0 for a flat Universe (required by 
inflationary models), k=l for a closed and k=-l for an open Universe. In 
this section we shall outline the main Cosmological Measurements and the 
pertinent quantities, of interest to us in these Lectures. For more details 
we refer the reader to the literature 12 . 
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2.1. Model Independent (Geometric) Considerations 

An important quantity, which we shall make extensive use of in the fol- 
lowing, when we use astrophysical data to constrain theoretical models, is 

the so-called Luminosity Distance, dL, defined as: = yj^p, where C is 
the energy per unit time emitted by the source, at the source's rest frame, 
and F is the flux measured by detector, i.e. the energy per unit time per 
unit area measured by the detector. Taking into account the time dilation 
induced by the expansion of the Universe: c>idetcctor = (<5i)sourcc(l + z), as 
well as as the cosmic red-shift, i.e. the fact that there is a reduced energy of 
photons at the detector as compared with that at emission from the source, 
we obtain: T = -. — ,. , , 2 . 

Another very commonly used quantity in Astrophysics is the Angular 
Diameter, which is defined as follows: A celestial object (cluster of galaxies 
etc.) has proper diameter D at r = r\ and emits light at t = t\. The 
observed angular diameter by a detector at t — to is: <5 = a J^y ri ■ From 
this one defines the Angular Diameter Distance: dA = -j = a{ti)r\ = 
d L (l + z) 2 . 

The Horizon Distance (beyond which light cannot reach us) is defined 
as: ds 2 = = dt 2 - a 2 (t) {j^pi +r 2 dQ 2 ^j . For radial motion of light, 
pertinent to most observations, along null geodesies ds 2 = 0, we have: 

So W) = So" 7& from Which dH ^ fl W So" V97r = a(t) Si $r y In 
Standard Cosmology dn ~ ^Ago < oo due to the finite age of the Universe, 
i.e. there is an Horizon. 

The above quantities are related among themselves 12 , as follows from 
the cosmic redshift phenomenon, the fact that photons follow null geodesies 
ds 2 = etc. These leads to relations among H ,dL and the redshift z, 
which are model independent, and they follow from pure geometrical con- 
siderations relying on the assumption of a RW homogeneous and isotropic 
cosmology. In the next subsection we discuss how a specific dynamical 
model of the Universe affects the cosmological measurements. In partic- 
ular, as we shall show, model dependence is hidden inside the details of 
the dependence of the Hubble parameter on the various components of the 
Universe's energy budget. This property is a consequence of the pertinent 
dynamical equations of motion of the gravitational field. 
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2.2. Cosmological Measurements: Model Dependence 

Within the standard General-Relativistic framework, according to which 
the dynamics of the gravitational field is described by the Einstein-Hilbcrt 
action, the gravitational (Einstein) equations in a Universe with cosmo- 
logical constant A read: R^v — %9 '^v R + g^uA- — SkGnThv , where Gat is 
(the four-dimensional) Newton's constant, Too = p is the energy density 
of matter, and Ta — a 2 (t)p with p = the pressure, and we assumed that 
the Universe and matter systems behave like ideal fluids in a co-moving 
cosmological frame, where all cosmological measurements are assumed to 
take place. From the RW metric (2.1), we arrive at the Friedman equation: 

' \ 2 k 

- + 3^ - A = 8nG N p (2.2) 
a J a z 



From this equation one obtains the expression for the Critical density (i.e. 
the total density required for flat k = A = Universe): p c = s Jq n (f ) • 

From the dynamical equation (2.2) one can obtain various relations be- 
tween the Hubble parameter H(z), the luminosity distance di, the decel- 
eration parameter q(z) and the energy densities p at various epochs of the 
Universe. For instance, for matter dominating flat (k — 0) Universes with 
A > and various (simple, ^-independent) equations of state p — Wip , 
(w r = 1/3 (radiation), w m = (matter-dust), w\ — — 1 (cosmological 
constant (de Sitter)) we have for the Hubble parameter: 

H(z) = H + zf {1+w ^ (2.3) 

with the notation: f2j = — , i = r(adiation), m(attcr), A, ... 

For the deceleration parameter we have at late eras, where radiation 

is negligible : q(z) = -ffi = (^j) (|^m(l + zf - fl A ) , with q = 
\£l m — Oa- Thus, it becomes evident that A acts as "repulsive" grav- 
ity, tending to accelerate the Universe currently, and eventually dominates, 
leading to an eternally accelerating de Sitter type Universe, with a future 
cosmic horizon. At present in the data there is also evidence for past de- 
celeration (q(z) > , for some z > z* > 0), which is to be expected if the 
dark energy is (almost) constant, due to matter dominance in earlier eras: 

q(z) > o (i + zf > 2n A /n m => z> z* = (|^) 1/3 - 1 . 

Finally, for the luminosity distance we obtain the important relation: 
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d L = (l + z)J t 



fto a(t )dt 
ti a(t) 




, from which 



(2.4) 



We shall use this relation in the following, in order to constrain various 
theoretical cosmological models by means of astrophysical observations. 

2.3. Supernovae la Measurements of Cosmic Acceleration 

Type la Supernovae (SNe) behave as Excellent Standard Candles, and thus 
can be used to measure directly the expansion rate of the Universe at high 
redshifts (z > 1) and compare it with the present rate, thereby provid- 
ing direct information on the Universe's acceleration. SNe type la are very 
bright objects, with absolute magnitude M ~ 19.5, typically comparable to 
the brightness of the entire host galaxy! This is why they can be detected 
at high redshifts z <~ 1, i.e. 3000 Mpc, lpc ~ 3x 10 16 m. Detailed studies of 
the luminosity profile 1,2 of each SNe suggests a strong relation between the 
width of the light curve and the absolute luminosity of SNe. This implies 
an accurate determination of its absolute luminosity. For each supernova 
one measures an effective (rest frame) magnitude in blue wavelength band, 
m e Jf , which is then compared with the theoretical expectation (depend- 
ing on the underlying model for the Universe) to yield information on the 
various fP. The larger the magnitude the dimmer the observed SNe. 

To understand the pertinent measurements recall the relation between 
the observed (on Earth) and emitted wavelengths A Q b s = (1 + z)X cm i t , 
as a result of the cosmic rcdshift phenomenon. In a magnitudc-redshift 
graph, if nothing slowed down matter blasted out of the Big Bang, we 
would expect a straight line. The data from High-redshift (z <~ 1) SNe 
la, showed that distant SNe lie slightly above the straight line. Thus they 
are moving away slower than expected. So at those early days (z ~ 1) 
the Universe was expanding at a slower rate than now. The Universe ac- 
celerates today! In such measurements, one needs the Hubble-Constant- 



Free Luminosity Distance: T>l(z; Dm, ^a) = ^f-^L, = y j§p, with L 
the Intrinsic Luminosity of source, F the Measured Flux. In Fried- 
man models T>l is parametrically known in terms of Qmj^a- An im- 
portant quantity used in measurements is the Distance Modulus m - 



Magnitude of the Source, M the Absolute Magnitude, and M. = M — 
51ogi? + 25 the fit parameter. Comparison of theoretical expectations 




M, where m = M + 25 + 51og 




= M + 51ogI?L, with m= Apparent 
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with data restricts JT2m, ^A- An important point to notice is that for fixed 
redshifts z the eqs. T>l(z;Qm^a) =constant yields degeneracy curves 
C in the fi-plane, of small curvature to which one associates a small 
slope, with the result that even very accurate data can at best select a 
narrow strip in £1- plane parallel to C. The results (2004) are summa- 
rized in figure 2 In the early works (1999) it was claimed that the best 
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Figure 2. Supernovae (and other) measurements on the Universe's energy budget. 



fit model, that of a FRW Universe with matter and cosmological con- 
stant for z < 3 (where the SNe data are valid) yields the following val- 
ues: 0.8D,m — 0.6Oa — —0.2 ± 0.1 , for < 1.5. Assuming a flat model 
(k=0) the data imply: Q^ at = 0.28±^ (la stat)t° Q a l (identified syst.), 
that is the Universe accelerates today 

q = \si M -n A ~ -0.6 < o 

Further support on these results comes, within the SNe measurement frame- 
work, from the recent (> 2004) discovery 2 , by Hubble Space Telescope, 
ESSENCE and SNLS Collaborations, of more than 100 high-z (2 > z > 1) 
supernovae, pointing towards the fact that for the past 9 billion years the 
energy budget of the Universe is dominated by an approximately constant 
dark energy component. 
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2.4. CMB Anisotropy Measurements by WMAP1,3 

After three years of running, WMAP provided a much more detailed pic- 
ture of the temperature fluctuations than its COBE predecessor, which 
can be analyzed to provide best fit models for cosmology, leading to severe 
constraints on the energy content of various model Universes, useful for par- 
ticle physics, and in particular supersymmetric searches. Theoretically 12 , 
the temperature fluctuations in the CMB radiation are attributed to: (i) 
our velocity w.r.t cosmic rest frame, (ii) gravitational potential fluctuations 
on the last scattering surface (Sachs- Wolf effect), (iii) Radiation field fluc- 
tuations on the last scattering surface, (iv) velocity of the last scattering 
surface, and (v) damping of anisotropics if Universe re-ionizes after decou- 
pling. A Gaussian model of fluctuations 12 , favored by inflation, is in very 
good agreement with the recent WMAP data (see figure 3). The perfect 
fit of the first few peaks to the data allows a precise determination of the 
total density of the Universe, which implies its spatial flatness. The various 

Angular scale (deg) 
90 2 0.5 0.2 
7000 E ' ' ' ' : 
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Multipole moment I 

Figure 3. Red points (larger errors) are previous measurements. Black points (smaller 
errors) are WMAP measurements (G. Hinshaw, et al. arXiv:astro-ph/0302217). 

peaks in the spectrum of fig. 3 contain interesting physical signatures: 

(i) The angular scale of the first peak determines the curvature (but not 
the topology) of the Universe. 

(ii) The second peak -truly the ratio of the odd peaks to the even peaks- 
determines the reduced baryon density. 

(iii) The third peak can be used to extract information about the 
dark matter (DM) density (this is a model-dependent result, though - 
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standard local Lorentz invariance assumed, see discussion in later sections 
on Lorentz- violating alternative to dark matter models). 

The measurements of the WMAP 6 on the cosmological parameters of 
interest to us here are given in 6 , and reviewed in 7 . The WMAP results 
constrain severely the equation of state p = wp (p = pressure), pointing 
towards w < —0.78, if one fits the data with the assumption —1 < w (we 
note for comparison that in the scenarios advocating the existence of a 
cosmological constant one has w = —1). Many quintessence models can 
easily satisfy the criterion — 1 < w < —0.78, especially the supersymmetric 
ones, which we shall comment upon later in the article. Thus, at present, 
the available data are not sufficient to distinguish the cosmological constant 
model from quintessence (or more generally from relaxation models of the 
vacuum energy). The results lead to the chart for the energy and matter 
content of our Universe depicted in figure 1, and are in perfect agreement 
with the Supernovae la Data 1 . The data of the WMAP satellite lead to 
a new determination of fitotai = 1.02 ± 0.02, where fitotai = Ptotai/pc, due 
to high precision measurements of secondary (two more) acoustic peaks as 
compared with previous CMB measurements (c.f. figure 3). Essentially 
the value of f2 is determined by the position of the first acoustic peak in 
a Gaussian model, whose reliability increases significantly by the discovery 
of secondary peaks and their excellent fit with the Gaussian model 6 . 

Finally we mention that the determination of the cosmological param- 
eters by the WMAP team 6 , after three years of running, favors, by means 
of best fit procedure, spatially flat inflationary models of the Universe 13 . 
In general, WMAP gave values for important inflationary parameters, such 
as the running spectral index, n s (k), of the primordial power spectrum of 
scalar density fluctuations 8% 14 P (k) = |<5j;| 2 . The running scalar spectral 
index n s (k) is n s (k) — , where k is the co-moving scale. Basically 

inflation implies n s = 1. WMAP measurements yield n s = 0.96, thus favor- 
ing Gaussian primordial fluctuations, as predicted by inflation. For more 
details we refer the reader to the literature 6 ' 7 . 

2.5. Baryon Acoustic Oscillations (BAO) 

Further evidence for the energy budget of the Universe is obtained by De- 
tection of the baryon acoustic peak in the large-scale correlation function 
of SDSS luminous red galaxies 4 . The underlying Physics of BAO can be 
understood as follows: Because the universe has a significant fraction of 
baryons, cosmological theory predicts that the acoustic oscillations (CMB) 
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in the plasma will also be imprinted onto the late-time power spectrum of 
the non-relativistic matter: from an initial point perturbation common to 
the dark matter and the baryons, the dark matter perturbation grows in 
place while the baryonic perturbation is carried outward in an expanding 
spherical wave. At recombination, this shell is roughly 150 Mpc in radius. 
Afterwards, the combined dark matter and baryon perturbation seeds the 
formation of large-scale structure. Because the central perturbation in the 
dark matter is dominant compared to the baryonic shell, the acoustic fea- 
ture is manifested as a small single spike in the correlation function at 150 
Mpc separation 4 . 

The acoustic signatures in the large-scale clustering of galaxies yield 
three more opportunities to test the cosmological paradigm with the early- 
universe acoustic phenomenon: 

(1) They would provide smoking-gun evidence for the theory of gravita- 
tional clustering, notably the idea that large-scale fluctuations grow by 
linear perturbation theory from z <~ 1000 to the present; 

(2) they would give another confirmation of the existence of dark matter 
at z ~ 1000, since a fully baryonic model produces an effect much larger 
than observed; 

(3) they would provide a characteristic and reasonably sharp length scale 
that can be measured at a wide range of redshifts, thereby determining 
purely by geometry the angular-diameter-distance-redshift relation and the 
evolution of the Hubble parameter. 

In the current status of affairs of the BAO measurements it seems that 
there is an underlying-thcoretical-model dependence of the interpretation 
of the results, as far as the predicted energy budget for the Universe is 
concerned. This stems from the fact that for small deviations from Q m = 
0.3, = 0.7, the change in the Hubble parameter at z = 0.35 is about half 
of that of the angular diameter distance. Eisenstein et al. in 4 modelled 
this by treating the dilation scale as the cubic root of the product of the 
radial dilation times the square of the transverse dilation. In other words, 
they defined 

,1/3 / \ 1 ' 2 



D v {z) = 



D M {zY 



H(z) 



H = H Q [J2^ + z) 3il+W ' ) ) (2-5) 



where H(z) is the Hubble parameter and Dm{z) is the co-moving angular 
diameter distance. As the typical redshift of the sample is z = 0.35, we 
quote the result 4 for the dilation scale as Dy(0.35) = 1370 ± 64Mpc. 
The BAO measurements from Large Galactic Surveys and their results 
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Figure 4. Left: Residual magnitude versus redshift for supernovae from the 'gold' and 
the SNLS datasets for various cosmological models. Right: The Hubble-parameter vs. 
redshift relation for these models and observational data. The bands represent 68% 
confidence intervals derived by the SN analysis for the standard ACDM, the super- 
horizon (no DE) and the Q-cosmology models. The black rectangle shows the WMAP3 
estimate for Ho, the squares show the measurements from SDSS galaxies, the triangles 
result from high- z red galaxies, and the circles correspond to a combined analysis of 
supernovae data (from [16]). 

for the dark sector of the Universe are consistent with the WMAP data, 
as far as the energy budget of the Universe is concerned, but the reader 
should bear in mind that they based their parametrization on standard 
FRW cosmologies, so the consistency should be interpreted within that 
theory framework. 

2.6. Measuring H(z): an important constraint on models 

The previous results, based on SNe, CMB and BAO measurements, relied 
on the standard FRW Cosmological model for the Universe as the underly- 
ing theory. However, in modern approaches to (quantum) gravity, such as 
brane and string theories, the underlying dynamics may no longer be de- 
scribed by the simple Einstein-Hilbert action. One may have extra fields, 
such as the dilaton or moduli fields in theories with extra dimensions, plus 
higher order curvature terms which could become important in the early 
Universe. Moreover, there have been suggestions in the literature 15 that 
the claimed Dark Energy may not be there, but simply be the result of 
temperature fluctuations in a (flat) Universe filled with matter Qm = 1 
( "super- horizon model"). All such alternative theories should be tested 
against each one of the above-mentioned categories of measurements to- 
gether with an independent measurement of the behavior of the Hubble 
parameter vs. the redshift H(z), the latter coming from large galactic 
surveys. This latter measurement provides an important constraint which 
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could differentiate among the potential Dark Energy (DE)/Dark Matter 
(DM) models and their alternatives. This extra measurement has the po- 
tential of ruling out alternative models (to DM and DE) that otherwise 
fit the supernova data alone (in a m cff vs z plot). This happens, for in- 
stance, with the super-horizon model of 15 . I mention in passing that other 
non-equilibrium stringy cosmologies 10 , with relaxing to zero dark energy 
(quintessence- like due to the dilaton field) survive at present this constraint, 
as illustrated in figure 4. For more details I refer the reader to 16 and ref- 
erences therein. 

2.7. Cosmic Coincidence and Cosmological Constant Issues 

There may be several possible explanations regarding the Dark Energy part 
of the Universe's energy budget: 

(i) The dark energy is an "Honest" Cosmological Constant A ~ lO -122 -/^, 
strictly unchanging through space and time. This has been the working 
hypothesis of many of the best fits so far, but I stress it is not the only 
explanation consistent with the data. 

(ii) Quintessence: The Cosmological constant is mimicked by a slowly- 
varying field, 0, whose time until it reaches its potential minimum is (much) 
longer than the Age of Universe. Simplest Quintessence models assume 
exponential potentials V(<p) <~ e^. In such a case the pertinent equation 

V(6) 

of state reads: w — — . For 6 = — 2 mi one has a rclaxing-to-zcro 

{ -^-+v{4>) 

vacuum energy A(t) ~ const/i 2 (in Planck units), of the right order of 
magnitude today. Such a situation could be met 10 in some models of 
string theory, where the role of the quintessence field could be played by 
the dilaton 17 , i.e. the scalar field of the string gravitational multiplct. 

(iii) Einstein-Friedman model is incorrect, and one could have modifica- 
tions in the gravitational law at galactic or supergalactic scales. Models of 
this kind have been proposed as alternatives to dark matter, for instance 
Modified Newtonian Dynamics (MOND) by Milgrom 18 , and its field theory 
version by Bekenstein 19 , known as Tensor- Vector- Scalar (TeVeS) theory, 
which however, is Lorentz Violating, as it involves a preferred frame. Other 
modifications from Einstein theory, which however maintain Lorentz invari- 
ancc of the four-dimensional world, could be brane models for the Universe, 
which are characterized by non-trivial, and in most cases time dependent, 
vacuum energy. It should be noted that such alternative models may lead 
to completely different energy budget 20 > 21 . We shall discuss one such case 
of a non-critical string inspired (non-equilibrium, relaxation) cosmology (Q- 
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cosmology) in a subsequent section, where we shall see that one may still 
fit the astrophysical data with exotic forms of "dark matter" , not scaling 
like dust with the redshift at late epochs, and different percentages of dark 
(dilaton quintessence) energy (c.f. also fig. 4). 

Given that from most of the standard best fits for the Universe it follows 
that the energy budget of our Cosmos today is characterized by 73 — 74% 
vacuum energy, i.e. an energy density of order p vac ~ (1CP 3 eU) 4 = 
1CP 8 erg/cm 3 , and about 27 — 26% matter (mostly dark), this implies 
the Coincidence Problem: "The vacuum energy density today is approxi- 
mately equal (in order of magnitude) to the current matter density." As the 
Universe expands, this relative balance is lost in models with a cosmologi- 
cal constant, such as the standard ACDM model, since the matter density 
scales with the scale factor as -U^ = — a a 3 . In this framework, at early 
times we have that the Vacuum Energy is much more suppressed as com- 
pared with that of Matter and Radiation, while at late times it dominates. 
There is only one brief epoch for which the transition from domination of 
one component to the other can be witnessed, and this epoch, according 
to the ACDM model, happened to be the present one! This calls for a 
microscopic Explanation, which is still lacking. 

The smallness of the value of the Dark Energy today is another big mys- 
tery of particle physics. For several years the particle physics community 
thought that the vacuum energy was exactly zero, and in fact they were 
trying to devise microscopic explanations for such a vanishing by means of 
some symmetry. One of the most appealing, but eventually failed in this 
respect, symmetry justifications for the vanishing of the vacuum energy 
was that of supersymmetry (SUSY): if unbroken, supersymmetry implies 
strictly a vanishing vacuum energy, as a result of the cancelation among 
boson and fermion vacuum-energy contributions, due to opposite signs in 
the respective quantum loops. However, this cannot be the correct ex- 
planation, given that SUSY, if it is to describe Nature, must be broken 
below some energy scale M susy , which should be higher than a few TeV, 
as partners have not been observed as yet. In broken SUSY theories, in 
four dimensional space times, there are contributions to vacuum energy 
Pvac oc^ ftM 4 usy ~ (few TeV) 4 , which is by far greater than the observed 
value today of the dark energy A - 1(T 122 Mp, M P ~ 10 19 GeV. Thus, 
SUSY does not solve the Cosmological Constant Problem, which at present 
remains one of the greatest mysteries in Physics. 

In my opinion, the smallness today of the value of the "vacuum" energy 
density might point towards a relaxation problem. Our world may have 
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not yet reached equilibrium, from which it departed during an early-epoch 
cosmically catastrophic event, such as a Big Bang, or — in the modern ver- 
sion of string/brane theory — a collision between two brane worlds. This 
non equilibrium situation might be expressed today by a quintessence like 
exponential potential exp (</>) , where </> could be the dilaton field, which in 
some models 10 behave at late cosmic times as <j> ~ — 21nt. This would 
predict a vacuum energy today of order 1/i 2 , which has the right order 
of magnitude, if t is of order of the Age of the Universe, i.e. t ~ 10 60 
Planck times. Supersymmetry in such a picture may indeed be a symmetry 
of the vacuum, reached asymptotically, hence the asymptotic vanishing of 
the dark energy. SUSY breaking may not be a spontaneous breaking but 
an obstruction, in the sense that only the excitation particle spectrum has 
mass differences between fermions and bosons. To achieve phenomenologi- 
cally realistic situations, one may exploit 22 the string/brane framework, by 
compactifying the extra dimensions into manifolds with non-trivial "fluxes" 
(these are not gauge fields associated with electromagnetic interactions, but 
pertain to extra-dimensional unbroken gauge symmetries characterizing the 
string models). In such cases, fermions and bosons couple differently, due 
to their spin, to these flux gauge fields (a sort of generalized "Zeeman" ef- 
fects). Thus, they exhibit mass splittings proportional to the square of the 
"magnetic field", which could then be tuned to yield phenomenologically 
acceptable SUSY-splittings, while the relaxation dark energy has the cos- 
mologically observed small value today. In such a picture, SUSY is needed 
for stability of the vacuum, although today, in view of the landscape scenar- 
ios for string theory, one might not even have supersymmetric vacua at all. 
However, there may be another reason why SUSY could play an important 
physical role, that of dark matter. I now come to discuss this important 
issue, mainly from a particle physics perspective. 

3. Dark Matter (DM) 

In this section I will discuss issues pertaining to dark matter and super- 
symmetry. I will first make the case for Dark Matter, starting historically 
from discrepancies concerning rotational curves of galaxies. Then I will 
move to describe possible candidates, and based on standard models for 
cosmology to exclude many of them, by means of WMAP data, arguing 
that supersymmetric dark matter remains compatible with such data. I 
will again emphasize, however, the model dependence of such conclusions. 
Then I will proceed to discuss supersymmetric particle physics constraints 
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in various frameworks by describing the underlying general framework for 
calculating thermal dark matter relics and compare them with WMAP 
data. For a more complete discussion on direct searches for dark matter 
the reader is referred to 23 , and references therein. 

3.1. The Case for DM 

Dark Matter (DM) is defined as a Non luminous massive matter, of un- 
known composition, that does not emit or reflect enough electromagnetic 
radiation to be observed directly, but whose presence can be inferred from 
gravitational effects on visible matter. Observed phenomena consistent 
with the existence of dark matter are: 

(i) rotational speeds of galaxies and orbital velocities of galactic clusters, 

(ii) gravitational lensing of background objects by galaxy clusters such as 
the Bullet cluster of galaxies, and 

(iii) the temperature distribution of hot gas in galaxies and clusters of 
galaxies. 

(iv) As we have seen, DM also plays a central role in structure formation 
and galaxy evolution, and has measurable effects on the anisotropy of the 
cosmic microwave background, especially the third peak in the anisotropy 
spectrum (c.f. fig. 3). 

400 
„ 350 



J 250 
b 

8 200 
> 150 
100 
50 

5 10 15 20 25 30 

R kpc 

Figure 5. Collage of Rotational Curves of nearby spiral galaxies obtained by combining 
Doppler data from CO molecular lines for the central regions, optical lines for the disks, 
and HI 21 cm line for the outer (gas) disks. Graph from Y. Sophue and V. Rubin (Annual 
Review of Astronomy and Astrophysics, Volume 31 (c)2001, 127). 

Historically, the first evidence for DM came 24 from discrepancies con- 
cerning the Rotational Curves (RC) of Galaxies. If all matter were luminous 
then the rotational speed of the galactic disc would fall with the (radial) 
distance r from the center as v(r) <~ r -1 / 2 but observations show that 
v(r) <~ const, as seen clearly in figure 5, where the rotation velocity in units 




16 



of km s _1 is plotted vs galactocentric radius R in kiloparsecs (kpc); 1 kpc 
« 3000 light years. It is seen that the RCs are flat to well beyond the 
edges of the optical disks (~ 10 kpc). Further Evidence for DM is provided 
by the Matter oscillation spectrum in galaxies, depicted in figure 6. The 
observed spectrum does not have the pronounced wiggles predicted by a 
baryon-only model, but it also has significantly higher power than does the 
model. In fact, A 2 = k 3 P(k)/(2ir 2 ) , which is a dimensionless measure of 
the clumping, never rises above one in a baryon-only model, so we could 
not see any large structures (clusters, galaxies, people, etc.) in the universe 
in such a model 25 . 
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Figure 6. Power spectrum of matter fluctuations (red curve, with wiggles) in a theory 
without dark matter as compared to observations of the galaxy power spectrum. 

However, at this stage we should mention the alternatives to Dark Mat- 
ter models, the MOND 18 , and its Lorentz- violating TeVeS field theory 
version 19 , which could also reproduce the rotational curves of galaxies, 
by assuming modified Newtonian dynamics at galactic scales for small 
gravitational accelerations, smaller than a universal value 7 < 70 ~ 
(200km sec _1 ) 2 /(10 kpc). MOND theories have been claimed to fit most 
of the rotational curves of galaxies (fig. 5), with few notable exceptions, 
though, e.g. the bullet cluster. It should be mentioned that TeVeS models, 
due to their preferred-cosmic- frame features, are characterized by "Aether" - 
Lorentz violating isotropic vector fields = (f(t), 0, 0, 0), A^A 11 = — 1, 
whose cosmic instabilities are also claimed 21 to reproduce the enhanced 
growth of perturbations observed in galaxies (c.f. fig. 6). In these lectures 
I will not discuss such models. It should be noted at this point that such 
issues, namely whether there arc dark matter particles or not, could be 
resolved in principle by particle physics searches at colliders or direct dark 
matter searches, which I will now come to. 
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3.2. Types of DM and Candidates 

From nucleosynthesis constraints we can estimate today the baryonic energy 
density contribution to be of order: Vt\, aryons = 0.045 ± .01, and this in fact 
is the dominant form of ordinary matter in the Universe. Thus, barring 
alternatives, 90% of the alleged matter content of the Universe seems to be 
dominated by DM of unknown composition at present. There are several 
dark matter candidates, which can be classified into two large categories 
depending on their origin and properties: 

(I) Astrophysical: (i) MAssive Compact Halo ObectS (MACHOS): Dwarf 
stars and Planets (Baryonic Dark Matter) and Black Holes, (ii) Non- 
luminous Gas Clouds. 

(II) Particles (Non-Baryonic Dark Matter): Weakly Interacting Massive 
Particles (WIMP), which might be the best candidates for DM: should not 
have electromagnetic or strong interactions. May have weak and gravita- 
tional interactions. WIMPS might include axions, neutrinos stable super- 
symmetric partners etc. If these WIMPS are thermal relics from the Big 
Bang then we can calculate their relic abundance today and compare with 
CMB and other astrophysical data. Non-thermal relics may also exist in 
some cosmological models but will not be the subject of our discussion in 
these lectures. 

There is an alternative classification of DM, depending on the energetics 
of the constituting particles: 

(i) Hot Dark Matter (HDM): form of dark matter which consists of particles 
that travel with ultra-relativistic velocities: e.g. neutrinos. 

(ii) Cold Dark Matter ( CDM): form of dark matter consisting of slowly 
moving particles, hence cold, e.g. WIMPS (stable supersymmetric particles 
(e.g. neutralinos etc.) or MACHOS. 

(iii) Warm Dark Matter (WDM): form of dark matter with properties 
between those of HDM and CDM. Examples include sterile neutrinos, light 
gravitinos in supergravity theories etc. 

Particle physics and/or astrophysics should provide candidates for DM 
and also explain the relic densities of the right order as that predicted by 
the data. Currently, the most favorite SUSY candidate for non baryonic 
CDM are neutralinos 26 \. These particles could be a WIMP if they are sta- 
ble, which is the case in models where they are the Lightest SUSY Particles 
(LSP) (with typical masses > 35 GeV). Most of supersymmetric model 
constraints come from the requirement that a neutralino is the dominant 
astrophysical DM, whose relic abundance can explain the missing Universe 
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mass problem. I mention at this stage that direct searches for \ involve, 
among others, the recoil of nucleons during their interaction with % in cryo- 
genic materials. In these lectures we shall concentrate mainly on colliders 
DM searches. I refer the reader to ref. 23 for direct DM searches and other 
pertinent terrestrial and extraterrestrial experiments. 

3.3. WIMP DM: thermal properties and relic densities 

In all the searches we shall deal with in the present work, which are also the 
most commonly studied in the literature, one makes the standard assump- 
tion that the dark matter particle, \, is a thermal relic of the Big Bang: 
when the early Universe was dense and hot, with temperature T m x , 
X was in thermal equilibrium; annihilation of \ an d its antiparticle \ m to 
lighter particles, XX — > ^ an d the inverse process 11 —>■ XX proceeded with 
equal rates 12 . As the Universe expanded and cooled down to a temperature 
T < m x , the number density of x dropped exponentially, n x <~ e~ m */ T . 
Eventually, the temperature became too low for the annihilation to keep up 
with the expansion rate and the species x 'froze out' with the cosmological 
abundance ( "relic" ) observed today. 

The time evolution of the number density n x (t) is determined by the 
Boltzmann equation 12 , 

dn x /dt + 3Hn x = ~{a A v) [{n x f - (n^) 2 ] , (3.6) 

where H is the Hubble expansion rate, n^ q the equilibrium number density 
and (<jav) is the thermally averaged annihilation cross section summed over 
all contributing channels. It turns out that the relic abundance today is 
inversely proportional to the thermally averaged annihilation cross section, 
il x h 2 ~ 1/{<jav). The situation is depicted in fig. 7. When the properties 
and interactions of the WIMP are known, its thermal relic abundance can 
hence be computed from particle physics' principles and compared with 
cosmological data. 

3.4. Hot and Warm DM Excluded by WMAP 

The WMAP / CMB results on the cosmological parameters discussed previ- 
ously disfavor strongly Hot Dark Matter (neutrinos), as a result of the new 
determination of the upper bound on neutrino masses. The contribution of 
neutrinos to the energy density of the Universe depends upon the sum of 
the mass of the light neutrino species 12,6 : 

^ - SB <"> 
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Figure 7. The full line is the equilibrium abundance; the dashed lines are the actual 
abundance after freeze-out. As the annihilation cross section (<ja_v) is increased, the 
WIMP stays in equilibrium longer, leading to a smaller relic density (from [12]). 

where the sum includes neutrino species that are light enough to decouple 
while still relativistic. 

The combined results from WMAP and other experiments 6 on the 
cumulative likelihood of data as a function of the energy density in neutrinos 
lead to 0„/i 2 < 0.0067 (at 95% confidence limit). Adding the Lyman a data, 
the limit weakens slightly 6 : Vl u h 2 < 0.0076 or equivalently (from (3.7)): 
J2i m vi < 0.69 eV, where, we repeat again, the sum includes light species of 
neutrinos. This may then imply an average upper limit on electron neutrino 
mass < m v > e < 0.23 eV. These upper bounds strongly disfavors Hot Dark 
Matter scenarios. 

Caution should be exercised, however, when interpreting the above 
WMAP result. There is the undcrlying-theoretical-model dependence of 
these results, which stems from the assumption of an Einstein-FRW Cos- 
mology, characterized by local Lorentz invariance. If Lorentz symmetry is 
violated, as, for instance, is the case of the TeVeS models alternative to 
DM, then neutrinos with (rest) masses of up to 2 eV could have an abun- 
dance of Q, u ~ 0.15 in order to reproduce the peaks in the observed CMB 
spectrum (fig. 3) 20 and thus being phenomenologically acceptable, at least 
from the CMB measurements viewpoint. 

At this juncture we note that another important result of WMAP is 
the evidence for early re-ionization of the Universe at redshifts z ~ 20. If 
one assumes that structure formation is responsible for re-ionization, then 
such early re-ionization periods are compatible only for high values of the 
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masses mx of Warm Dark Matter . Specifically, one can exclude models 
with mx < 10 KeV based on numerical simulations of structure formation 
for such models 27 . Such simulations imply that dominant structure forma- 
tion responsible for re-ionization, for Warm Dark Matter candidates with 
m x < 10 KeV, occurs at much smaller z than those observed by WMAP. 
In view of this, one can therefore exclude popular particle physics mod- 
els employing light gravitinos (m x — 0.5 KeV) as the Warm Dark Matter 
candidate. It should be noted at this stage that such structure formation 
arguments can only place a lower bound on the mass of the Warm Dark 
Matter candidate. The reader should bear in mind that Warm Dark Mat- 
ter with masses mx > 100 KeV becomes indistinguishable from Cold Dark 
Matter, as far as structure formation is concerned. 

3.5. Cold DM in Super symmetric Models: Neutralino 

After the exclusion of Hot and Warm Dark Matter, the only type of Dark 
matter that remains consistent with the recent WMAP results 6 is the Cold 
Dark Matter , which in general may consist of axions, superheavy particles 
(with masses ~ 10 14±5 GeV) 28,29 and stable supersymmetric partners. 
Indeed, one of the major and rather unexpected predictions of Supersym- 
metry (SUSY), broken at low energies Msusy ~ 0{1 TeV), while i?-parity 
is conserved, is the existence of a stable, neutral particle, the lightest neu- 
tralino (x), referred to as the lightest supersymmetric particle (LSP) 26 . 
Such particle is an ideal candidate for the Cold Dark Matter in the Uni- 
verse 26 . Such a prediction fits well with the fact that SUSY is not only 
indispensable in constructing consistent string theories, but it also seems 
unavoidable at low energies (~ 1 TeV) if the gauge hierarchy problem is to 
be resolved. Such a resolution provides a measure of the SUSY breaking 
scale M SUS Y ~ 0(1 TeV). 

This type of Cold Dark Matter will be our focus from now on, in 
association with the recent results from WMAP3 on relic densities 3 ' 6 . 
The WMAP3 results, combined with other existing data, yield for the 
baryon and matter densities (including dark matter) at 2<r-level: Vt m h 2 — 
0.1268±gS (matter) , 100f2 6 /i 2 = 2.233±g;g^ (baryons) . One assumes 
that CDM is given by the difference of these two. As mentioned already, 
in supersymmetric (SUSY) theories the favorite candidate for CDM is the 
lightest of the Neutralinos \ (SUSY CDM), which is stable, as being the 
Lightest SUSY particle (LSP) (There are cases where the stau or the sncu- 
trino can be the lightest supersymmetric particles. These cases are not fa- 
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vored 30 and hence are not considered). From the WMAP3 results 3 , then, 
assuming Sic dm — we can infer stringent limits for the neutralino \ 
relic density: 

0.0950 < Q x h 2 < 0.1117 , (3.8) 

It is important to notice that in this inequality only the upper limit is 
rigorous. The lower Limit is optional, given that there might (and probably 
do) exist other contributions to the overall (dark) matter density. It is 
imperative to notice that all the constraints we shall discuss in this review 
are highly model dependent. The results on the minimal SUSY extensions 
of the standard model 7 , for instance, cannot apply to other models such as 
superstring-inspired ones, including non equilibrium cosmologies, which we 
shall also discuss here. However, formally at least, most of the analysis can 
be extrapolated to such models, with possibly different results, provided the 
SUSY dark matter in such models is thermal. Before moving into such a 
discussion we consider it as instructive to describe briefly various important 
properties of the Neutralino DM. 

The Neutralino is a superposition of SUSY partner states. Its mass 
matrix in bino-wino-higgsinos basis tp® = (—iX 1 , —i\ 3 , ^h 2 ) ^ s gi yen 

by ' 



M 



N 



/Mi -m z s w cp mzsy/sp \ 

M 2 m z c w c f j -m z cwSf] 

-mz-swcp mzcwcp — ji 

\ m z s w sp -mzCwS/3 -/i / 



where Mi, M 2 : the U(l) and SU(2) gaugino masses, [i: higgsino mass 
parameter, s\y — sin8w, cw — cos 9\y, S/3 = sin/3, cp — cos/3 and tan/3 = 
v 2 /vi (vi.2 v.e.v. of Higgs fields #1,2)- The mass matrix is diagonalized by a 
unitary mixing matrix N, N*MnN^ — diag(m^o, m^o, m^o, m^o) , where 
to^o, i — 1, ...,4, are the (non- negative) masses of the physical neutralino 
states with m^o < ... < m^a. The lightest neutralino is then: 

xl = N^B + N 12 W + /V 13 i?i + N U H 2 . 

To calculate relic densities it is assumed that the initial number density 
of neutralinos \ particle in the Early Universe has been in thermal equilib- 
rium: interactions creating x usually happen as frequently as the reverse 
interactions which destroy them. Eventually the temperature of the ex- 



panding Universe drops to the order of the neutralino (rest) mass T ~ m 
In such a situation, most particles no longer have sufficient energy to create 
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neutralinos. Now neutralinos can only annihilate, until their rate becomes 
smaller than the Hubble expansion rate, H > r ann . Then, Neutralinos 
are being separated apart from each other too quickly to maintain equi- 
librium, and thus they reach their freeze-out temperature, Tf — m x /20, 
which characterizes this type of cold dark matter. 
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Figure 8. Basic Neutralino annihilations including stau co-annihilations in MSSM (from 
S. Kraml, Pramana 67, 597 (2006) [hep-ph/0607270]). 



In most neutralino relic density calculations, the only interaction cross 
sections that need to be calculated are annihilations of the type XX ~ > -X 
where \ 1S the lightest neutralino and X is any final state involving only 
Standard Model particles. However, there are scenarios in which other 
particles in the thermal bath have important effects on the evolution of the 
neutralino relic density. Such a particle annihilates with the neutralino into 
Standard Model particles and is called a co-annihilator (c.f. figure 8). In 
order for a particle to be an Effective co-annihilator, it must have direct 
interactions with the neutralino and must be nearly degenerate in mass: 
Such degeneracy happens in the Minimal Supersymmetric Standard Model 
(MSSM), for instance, with possible co-annihilators being the lightest stau, 
the lightest stop, the second-to-lightest neutralino or the lightest chargino. 
When this degeneracy occurs, the neutralino and all relevant co-annihilators 
form a coupled system. 

Without co-annihilations the evolution of a relic particle number den- 
sity, n, is governed, as mentioned previously, by a single-species Boltzmann 
equation (3.6). It should be noted that the relic-particle number density is 
modified by the Hubble expansion and by direct and inverse annihilations 
of the relic particle. The Relic particle is assumed stable, so relic decay is 
neglected. Also commonly assumed is time-reversal (T) invariance, which 
relates annihilation and inverse annihilation processes. In the presence of 
co-annihilators the Boltzmann equation gets more complicated but it can 
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(using n = Yh=i n r 



be simplified using stability properties of relic particle and co-annihilators 

=i«i): 

dn N 
— = - 3Hn - (<7H v ij) { n * n 3 - n T n T) ■ ( 3J ) 

To a very good approximation, one can use an effective single species Boltz- 
mann equation for this case if (crv) — j{&ijVij)-^-^. 

The Boltzmann equation (3.9) can be solved numerically, but in most 
cases even analytically. Details on how to solve the Boltzmann equation 
are given abundantly in the cosmology literature 12 and will not be re- 
peated here. We shall only outline the most important results that will 
be essential for our discussion in these lectures. One should determine the 
freeze-out temperature xp — m x /Tp: xp = In ^ 0,038 9 J^' J^* , with 
m p i the Planck mass, g the total number of degrees of freedom of the % 
particle (spin, color, etc.), the total number of effective relativistic de- 
grees of freedom at freeze-out, and the thermally averaged cross section 
is evaluated at the freeze-out temperature. For most CDM candidates, 
xp ~ 20. The total (co) annihilation depletion of neutralino number den- 
sity is calculated by integrating the thermally averaged cross section from 
freeze-out to the present temperature: 

h 2 _ f^]f_Jo 1 1-07 x K^GcY -1 

J(x F ) = / (crv)x- 2 dx . (3.10) 

where sq is the entropy density, g*s denotes the number of effective rel- 
ativistic d.o.f. contributing to the (constant) entropy of the universe and 
h is the reduced Hubble parameter: Ho = 100 h kmsec -1 Mpc -1 . This is 
the expression one compares with the experimental determination of the 
DM abundance via, e.g., WMAP data. It should be noted at this stage 
that the theoretical assumptions leading to the above results may not hold 
in general for all DM models and candidates: the missing non-baryonic 
matter in the universe may only partially, or not at all, consist of relic 
neutralinos. Also, as we shall discuss later on the article, in some off-shell, 
non-equilibrium relaxation stringy models of dark energy, the Boltzmann 
equation gets modified by off-shell, non-equilibrium terms as well as time- 
dependent dilaton-source terms. This leads to important modifications on 
the associated particle-physics models constraints. 
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Figure 9. Values of the quantity a a n allowed at 2<x level as a function of the DM mass. 

3.6. Model- Independent DM Searches in Colliders 

As we have discussed above, if dark matter comes from a thermal relic, 
its density is determined, to a large extent, by the dark matter annihila- 
tion cross section: cr(xx ~~ * SM SM). Indeed, as already mentioned, the 
present-day dark matter abundance is roughly inversely proportional to 
the thermally averaged annihilation cross section times velocity, fl x h 2 oc 
1 / (av) . This latter quantity can be conveniently expanded in powers of the 
relative dark matter particle velocity: 

= 5>&V aj > . ( 3 - n ) 

j 

Usually, only the lowest order non-negligible power of v dominates. For 
J = 0, such dark matter particles are called s-annihilators, and for J = 1, 
they are called p-annihilators; powers of J larger than 1 are rarely needed.: 
Figure 9 shows the constraint on the annihilation cross section as a 
function of dark matter mass that results from Eq. (3.8) 31 . The lower 
(upper) band of fig. 9 is for models where s-wave (p-wave) annihilation 
dominates. It is important to notice 31 that the total annihilation cross 
section cr an is virtually insensitive to dark-matter mass. This latter effect is 
due to the changing number of degrees of freedom at the time of freeze-out 
as the dark matter mass changes. It also points to cross sections expected 
from weak-scale interactions (around 0.8 pb for s-annihilators and 6pb for 
p-annihilators), hence implying the possibility that DM is connected to an 
explanation for the weak scale and thus WIMPs 31 . Such WIMPs exist 
not only in supcrsymmctric theories, of course, but in a plethora of other 
models such as theories involving extra dimensions and 'little Higgs' mod- 
els. The LHC and the ILC are specifically designed to probe the origin 
of the weak scale, so dark matter searches and future collider physics ap- 
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Figure 10. Left panel: Comparison between the photon spectra from the process 
e+e - — * 2\i + 7 in the explicit supersymmetric models defined in A. Birkedal, 
K. Matchev and M. Perelstein, Phys. Rev. D 70, 077701 (2004) (red/dark-gray) and 
the predicted spectra for a p-annihilator of the corresponding mass and K e (green/light- 
gray). Right panel: The reach of a 500 GeV unpolarized electron-positron collider with an 
integrated luminosity of 500 fb" 1 for the discovery of p-annihilator WIMPs, as a function 
of the WIMP mass M x and the e+e - annihilation fraction ft e . The 3 a (black) contour is 
shown, along with an indication of values one might expect from supersymmetric models 
(red dashed line, labelled 'SUSY'). Only statistical uncertainty is included. 



pear to be closely related. The next question one could ask is whether the 
above cross section could be turned, within a WIMP working hypothesis 
framework, into a model-independent signature at colliders. This ques- 
tion was answered in the affirmative in 31 . One introduces the parameter 
K e = a(xx ~ * e+e ~)/ cr (xX — * SM\SM) which relates dark matter annihi- 
lation processes to cross sections involving e + e _ in the final state. Using 
crossing symmetries to relate a(xx ~ * e + e~) to er(e + e - — > xx) an d co- 
linear factorization one can relate <r(e + e~ — * xx) to c(e + e~ — > xXl)i thus 
connecting astrophysical data on er an to the process e + e~ — > XXI- The 
resulting differential cross section reads 31 

-(e+e ^2 x + 7)- 



dxdcosO 

a K e a ari 1 + (I — x) 2 I 
16tt x sin 2 6> 



-2 Jo (2S' x + l) 2 1 




with a the appropriate fine structure constant, x = 2E 1 /^fs 1 angle be- 
tween photon and incoming electron, S x spin of WIMP, J is the dominant 
value of J in the velocity expansion of (3.11) (as discussed above, commonly 
J = dominates, s-annhilator DM). The accuracy of the method and its 
predictions are illustrated in fig. 10, where the left panel illustrates the 
results obtained using the formula (3.12), which are then compared with 
those of an exact calculation, based on a supersymmetric MSSM model, 
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with WIMP masses 225 GeV, whilst the right panel shows the expected 
reach in K e for a 500 GeV linear e + e~ collider as a function of the WIMP 
mass. As we observe from such comparisons the results of the method and 
of the exact calculation are in pretty good agreement. 

We note at this stage, however, that, although model independent, the 
above process is rarely the dominant collider signature of new physics within 
a given model. It therefore makes sense to look for model dependent pro- 
cesses at colliders, which we now turn to. In this last respect, it is important 
to realize 31 that a calculation of slepton masses is essential for comput- 
ing accurately relic abundances in theoretical models; without a collider 
measurement of the slepton mass, there may be a significant uncertainty 
in the relic abundance calculation. This uncertainty results because the 
slepton mass should then be allowed to vary within the whole experimen- 
tally allowed range. We mention here that measuring slepton masses at 
LHC is challenging due to W + W~ and tt production. However, as shown 
in 31 , it is possible through the study of di-lepton mass distribution m« 
in the decay channel X2 ~^ ^ ^Xi an d a l so a t the International Linear 
Collider (ILC). The reader is referred to the literature 31 for further details 
on these important issues. We are now ready to start our discussion on 
model-dependent DM signatures at LHC and future colliders. 

4. Model-Dependent WMAP SUSY Constraints 

We shall concentrate on DM signatures at colliders, using WMAP 1,3 data. 
To illustrate the underlying-theoretical-model dependence of the results we 
chose three representative theoretical models: (i) the mSUGRA (or con- 
strained MSSM model) 8 ' 7 , (ii) a heterotic string model with orbifold com- 
pactification 9 , and (iii) a non-critical (non-equilibrium) stringy cosmology 
(Q-cosmology) with running dilatons (implying a dilaton quintessence re- 
laxation model for dark energy at late eras) and off-shell terms 10 . 

4.1. Constrained MSSM/mS UGRA Model 

MSSM has too many parameters to be constrained effectively by data. To 
minimize the number of parameters one can "embed" this model by taking 
into account the gravity sector, which from a cosmological point of view is 
a physical necessity. Such an embedding in principle affects the dark en- 
ergy sector of the cosmology, and in fact the minimal Supergravity model 
(mSUGRA) 8 , used to yield the Constrained MSSM (CMSSM) predicts too 
large values of the cosmological constant at a quantum level, and hence it 
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should not be viewed as the physical model. Nevertheless, as far as DM 
searches are concerned, such models give a pretty good idea of how astro- 
physical data can be used to constrain particle physics models, and this is 
the point of view we take in this work. mSUGRA is the best studied model 
so far as far as constraints on supersymmetric models using astrophysical 
CMB data are concerned. A relatively recent review on such approaches is 
given in 7 , where we refer the reader for details and further material and 
references. In our presentation here we shall be very brief and concentrate 
only on the basic conclusions of such analyses. 

4.1.1. Basic Features: geometry of the parameter space 

Before embarking into a detailed analysis of the constraints of the minimal 
supersymmetric standard model embedded in a minimal supergravity model 
(CMSSM) 8 , we consider it useful to outline the basic features of these 
models, which will be used in this review. The embedding of SUSY models 
into the minimal supergravity (mSUGRA) model implies that there are five 
independent parameters: Three of them, the scalar and gaugino masses 
mo,mi/2 as well as the trilinear soft coupling Aq =, at the unification 
scale, set the size of the Supersymmetry breaking scale. In addition one 
can consider as input parameter tan/3 = , the ratio of the v.e.v's of 

the Higgses Hi and H\ giving masses to up and down quarks respectively. 
The sign ( signature) of the Higgsino mixing parameter fi is also an input 
but not its size which is determined from the Higgs potential minimization 
condition 7 . The parameter space of mSUGRA can be effectively described 
in terms of two branches: 

(i) An Ellipsoidal Branch (EB) of Radiative Symmetry Breaking, which 
exists for small to moderate values of tan/3 <, 7, where the loop correc- 
tions are typically small. One finds that the radiative symmetry breaking 
constraint demands that the allowed set of soft parameters ttlq and a combi- 
nation 7 m' 12 = f{mi/2, Aq, tan/3) lie, for a given value of (x, on the surface 
of an Ellipsoid. This places upper bounds on the sparticle masses for a 
given value of <f> = [i 2 jM 2 z + 1/4. 

(ii) Hyperbolic Branch (HB) of Radiative Symmetry Breaking. This 
branch is realized 32 for large values of tan/3 ^ 7, where the loop cor- 
rections to fi are significant. In this branch, (mo,m'^ 2 ) lie now on the 

surface of a hyperboloid: a ^Q ) ~ p™(Q B ) — Qo — & fixed value of 

the running scale, a, f3 constant functions of Mz, Aq. For fixed Aq, the 
mo, rax/2 ue on a hyperbola, hence they can get large for fixed \x or <E>. What 
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is interesting in the HB case is the fact that mo and/or m 1 / 2 can become 
very large, while much smaller values for \x can occur. 

(iia) A subset of HB is the so-called high zone 32 . In this case electroweak 
symmetry breaking (EWSB) can occur in regions where m and m 1 / 2 can 
be in the several TeV range, with much smaller values for the parame- 
ter \x which however is much larger than Mz- This has important con- 
sequences for phenomenology, as we shall see. In this zone the lightest 
of the neutralinos, xi; is almost a Higgsino having a mass of order fi. 
This is called inversion phenomenon since the LSP is a Higgsino rather 
a Bino. The inversion phenomenon has dramatic effects on the nature of 
the particle spectrum and SUSY phenomenology in this HB. Indeed, as 
we discussed above, in mSUGRA one naturally has co-annihilation with 
the sleptons when the neutralino mass extends to masses beyond 150- 
200 GeV with processes of the type (c.f. fig. 8): x^r -> £ a J,£ a Z,£ a h, 
£ a R £ b R £ a £ b , and £ a R £ b * -> ^,77,7^, ZZ, W+W~, hh, where 1 is es- 
sentially a f. Remarkably the relic density constraints can be satisfied 
on the hyperbolic branch also by co-annihilation. However, on the HB 
the co-annihilation is of an entirely different nature as compared with the 
stau co-annihilations discussed previously: instead of a neutralino-stau co- 
annihilation, and stau - stau in the HB one has co-annihilation processes 
involving the second lightest neutralino and chargino states 33 , X? ~ xf> 
followed by X? — X21X1 ~ Xi jXi — X2 • Some of the dominant processes that 
contribute to the above co-annihilation processes are 33 : x\Xi iX^Xi ~ * 
Uidi, SiVi, AW + , ZW + ,W + h and X1X1, X1X2 ~~ > 1 didi , W + W~ . Since 
the mass difference between the states xt an d X? is the smallest the X1X1" 
co-annihilation dominates. In such cases, the masses too m 1 / 2 may be 
pushed beyond 10 TeV, so that squarks and sleptons can get masses up to 
several TeV, i.e. beyond detectability limits of immediate future accelera- 
tors such as LHC. 

(iib) Except the high zone where the inversion phenomenon takes place the 
HB includes the so called Focus Point (FP) region 34 , which is defined as 
a region in which some renormalization group (RG) trajectories intersect 
(FP region would be only a point, were it not for threshold effects which 
smear it out). We stress that the FP is not a fixed point of the RG. 
The FP region is a subset of the HB limited to relatively low values of 
TO4/2 and values of \i close to the electroweak scale, Mz, while m can 
be a few TeV but not as large as in the high zone due to the constraints 
imposed by the EWSB condition. The LSP neutralino in this region is a 
mixture of Bino and Higgsino and the Higgsino impurity allows for rapid 
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s-channel LSP annihilations, resulting to low neutralino relic densities at 
experimentally acceptable levels. This region is characterized by mo in 
the few TeV range, low values of mi 1/2 << m and rather small values 
of \x close to Mz- The LSP neutralino in this case is a mixture of Bino 
and Higgsino and its Higgsino impurity is adequate to give rize to rapid 
s-channel LSP annihilations so that the neutralino relic density is kept low 
at experimentally acceptable values. Since \i is small the lightest chargino 
may be lighter than 500 GeV and the FP region may be accessible to 
future TeV scale colliders. Also due to the relative smallness of mi/ 2 in 
this region gluino pair production may occur at a high rate making the FP 
region accessible at LHC energies. 

It should be pointed out that, although the HB may be viewed as fine 
tuned, nevertheless recent studies 35 , based on a % 2 analysis, have indicated 
that the WMAP data, when combined with data on b — > and g^ — 2, 
seem to favor the Focus Point HB region and the large tan (3 neutralino 
resonance annihilation of mSUGRA. 

4.1.2. Muon's anomaly and SUSY detection prospects 

Undoubtedly one of the most significant experimental results of the last 
years is the measurement of the anomalous magnetic moment of the 
muon 36 . Deviation of its measured value from the Standard Model (SM) 
predictions is evidence for new physics with Supersymmetry being the 
prominent candidate to play that role. Adopting Supersymmetry as the 
most natural extension of the SM, such deviations may be explained and 
impose at the same time severe constraints on the predictions of the avail- 
able SUSY models by putting upper bounds on sparticle masses. Therefore 
knowledge of the value of — 2 is of paramount importance for Supersym- 
metry and in particular for the fate of models including heavy sparticles 
in their mass spectrum, as for instance those belonging to the Hyperbolic 
Branch. 

Unfortunately the situation concerning the anomalous magnetic mo- 
ment is not clear as some theoretical uncertainties remain unsettled as yet. 
Until last year, as far as I am aware of, there were two theoretical estimates 
for the difference of the experimentally measured 36 value of = (g tl — 2)/2 
from the theoretically calculated one within the SM 37 , 

• Estimate (I) a° xp - af* = 1.7(14.2) x 10~ 10 [0.4(15.5) x 10~ 10 ] 

• Estimate (II) a^-a™ = 24.1(14.0) x lO" 10 [22.8(15.3) x lO" 10 ] 
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In (I) the r-decay data are used in conjunction with Current Algebra while 
in (II) the e~e + — > Hadrons data are used in order to extract the photon 
vacuum polarization which enters into the calculation of — 2. Within 
square brackets are updated values of Ref . 37 1 . Estimate (I) is considered 
less reliable since it carries additional systematic uncertainties and for this 
reason in many studies only the Estimate (II) is adopted. Estimate (II) in- 
cludes the contributions of additional scalar mesons not taken into account 
in previous calculations. 

In order to get an idea of how important the data on the muon anomaly 
might be we quote Ref. 32 where both estimates have been used. If Estimate 
(II) is used at a 1.5a range much of the HB and all of the inversion region 
can be eliminated. In that case the usually explored region of SUSY in 
the EB is the only one that survives, which, as we shall discuss below, can 
be severely constrained by means of the recent WMAP data. On the other 
hand, Estimate (I), essentially implies no difference from the SM value, and 
hence, if adopted, leaves the HB, and hence its high zone (inversion region), 
intact. In such a case, SUSY may not be detectable at colliders, at least in 
the context of the mSUGRA model, but may be detectable in some direct 
dark matter searches, to which we shall turn to later in the article. 

For the above reasons, it is therefore imperative to determine unam- 
biguously the muon anomalous magnetic moment g^ — 2 by reducing the 
errors in the leading order hadronic contribution, experimentally, and im- 
proving the theoretical computations within the standard model. In view 
of its importance for SUSY searches, it should also be necessary to have 
further experiments in the future, that could provide independent checks of 
the measured muon magnetic moment by the E821 experiment 36 . Quite 
recently (2006) a new measurement 11 of the g^ — 2 became available, which 
shows a clear discrepancy from the theoretically calculated Standard Model 
prediction by 3.4 a 

1.91 x 10~ 9 < Aa M < 3.59 x 10~ 9 (4.13) 

thereby pointing towards the elimination of the inversion region of the HB 
of mSUGRA, according to the above discussion. 



x Due to the rapid updates concerning — 2 the values of — a^ M used in previous 
works quoted in this article may differ from those appearing above. 
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4.1.3. WMAP mSUGRA Constraints in the EB 

After the first year of running of WMAP, there have been two independent 
groups working on this update of the CMSSM in light of the WMAP data, 
with similar results 38 ' 39 and below we summarize the results of 38 in fig. 11 
for some typical values of the parameters tan/? and signature of fi. In such 
analyses one plots tuq vs. mi/2, taking into account the calculated relic 
abundance of neutralinos in the model and constraining it by means of the 
WMAP results (3.8). Details are given in 7 and references therein. 



Ian B = 10 , 11 > tan B = 10 , n < 

! is. 300 ■ i n , r , 
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Figure 11. mSUGRA/CMSSM constraints after WMAP from Ref. [38]. The Dark 
Blue shaded region is favored by WMAP1 (0.094 < tl x h 2 < 0.129 ). Turquoise shaded 
regions have 0.1 < Q. x h 2 < 0.3. Brick red shaded regions are excluded because LSP is 
charged. Dark green regions are excluded by b — » sj. The Pink shaded region includes 
2 — a effects of g M — 2. Finally, the dash-dotted line represents the LEP constraint on 
e mass. 

For the LSP, the lightest of the charginos, stops, staus and Higgses 
the upper bounds on their masses of order of a few hundreds of GeV 7 , for 
various values of the parameter tan j3 , if the new WMAP determination 3,6 
of the Cold Dark Matter (3.8) and the 2a bound 149 < lCT 11 af USY < 
573 of E821 is respected. The lightest of the charginos has a mass whose 
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upper bound is w 550 GeV , and this is smaller than the upper bounds 
put on the masses of the lightest of the other charged sparticles, namely 
the stau and stop. Hence the prospects of discovering CMSSM at a e+e~ 
collider with center of mass energy y/s = 800 GeV, are not guaranteed. 
Thus, a center of mass energy of at least y/s w 1.1 TcV is required to 
discover SUSY through chargino pair production. Note that in the allowed 
regions the next to the lightest neutralino, x', has a mass very close to 
the lightest of the charginos and hence the process e+e~ — > xx', with %' 
subsequently decaying to x + l + l~ or x + 2jets, is kinematically allowed for 
such large tan f3, provided the energy is increased to at least y/s = 860 GeV. 
It should be noted however that this channel proceeds via the f-channel 
exchange of a selectron and it is suppressed due to the heaviness of the 
exchanged sfermion. Therefore only if the center of mass energy is increased 
to y/s = 1.1 TeV supersymmetry can be discovered in a e + e~ collider 
provided it is based on the Constrained scenario 39 . 

An important conclusion, therefore, which can be inferred by inspecting 
the figures 11 is that the constraints implied by a possible discrepancy of 

— 2 from the SM value, as seems to be supported by the 2006 data 11 
(4.13), ( al USY k, 15.0 x 10~ 10 ), when combined with the WMAP restric- 
tions on CDM (neutralino) relic densities (3.8), imply severe restrictions on 
the available parameter space of the EB and lower significantly the upper 
bounds on the allowed neutralino masses m^. 



4.1.4. WMAP mSUGRA Constraints in the HB 

Despite the above-mentioned good prospects of discovering minimal SUSY 
models at future colliders, if the EB is realized, however, things may not 
be that simple in Nature, x 2 studies 35 of mSUGRA in light of the recent 
WMAP data has indicated (c.f. figure 12) that the HB/focus point region of 
the model's parameter space seems to be favored along with the neutralino 
resonance annihilation region for \i > and large tan/3 values. The favored 
focus point region corresponds to moderate to large values of the Higgs 
parameter [i 2 , and large scalar masses to in the several TeV range. The 
situation in case the HB is included in the analysis is depicted in figure 
13 32 , where we plot the too — mi« graphs, as well as graphs of mo, TO1/2 
vs the neutralino LSP mass. The neutralino density is that of the WMAP 
data. 

We stress again that, in case the high zone (inversion) region of the HB 
is realized, then the detection prospects of SUSY at LHC are diminished 



33 



mSugm with lanfi = 10, ■%=(!, jl>0 mSugm with tanfi = 56, A Q = 0, jl > 




1000 2000 3000 4000 5000 6000 1000 2000 3000 4000 5000 6000 

m (GeV) m (GeV) 

— m h =114.1GeV ■ LEP2 excluded — m h =114.1GeV ■ LEP2 excluded 

GENIUS CDMSII CDMS GENIUS CDMSII CDMS 



Figure 12. WMAP data seem to favor (-j^j < 4/3) (green) the HB/focus point region 
(moderate to large values of fi, large mo scalar masses) for almost all tan/3 (Left), as 
well as s - channel Higgs resonance annihilation (Right) for fi > and large tan/3. 



significantly, in view of the fact that in such regions slepton masses may lie 
in the several TeV range (see figure 13). Fortunately, as already mentioned, 
last years's — 2 data 11 (4.13) seem to exclude this possibility. 
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Figure 13. mo — mi« graph, and mo and vs - m x graphs, including the HB of 

mSUGRA. Such regions are favored by the WMAP data. 
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4.1.5. Expected Reach of LHC and Tevatron 

In view of the above results, an updated reach of LHC in view of the recent 
WMAP and other constraints discussed above (see figure 14) has been 
performed in 40 , showing that a major part of the HB, but certainly not its 
high zone (which though seems to be excluded by means of the recent <? M — 2 
data (4.13)), can be accessible at LHC. The conclusion from this study 40 



mSugra with w/ifi = 10, A u = 0, )i > 




Figure 14. Left: The updated Reach in (too, in 1/2) parameter plane of mSUGRA as- 
suming 100 fb~ 1 integrated luminosity. Red (magenta) regions are excluded by theoreti- 
cal (experimental) constraints. Right: Contours (in view of the uncertainties) of several 
low energy observables : CDM relic density (green), contour of = 114.1 GeV (red), 
contours o/a M 10 10 (blue) and contours of b — > s-y BF ('xlO 4 ) (magenta). 



is that for an integrated luminosity of 100 values of TO1/2 ~ 1400 GeV 
can be probed for small scalar masses mo, corresponding to gluino masses 
rrig ~ 3 TeV. For large mo, in the hyperbolic branch/focus point region, 
TO i/2 ~ 700 GeV can be probed, corresponding to m g ~ 1800 GeV. It is also 
concluded that the LHC (CERN) can probe the entire stau co-annihilation 
region and most of the heavy Higgs annihilation funnel allowed by WMAP 
data, except for some range of mo, m 1 / 2 in the case tan/3 <; 50. A similar 
updated reach study in light of the new WMAP data has also been done 
for the Tevatron 41 , extending previous analyses to large mo masses up to 
3.5 TeV, in order to probe the HB/focus region favored by the WMAP 
data 35 . Such studies (c.f. figure 15) indicate that for a 5cr (3tr) signal with 
10 (25) fb^ 1 of integrated luminosity, the Tevatron reach in the trilepton 
channel extends up to mi/ 2 ~ 190 (270) GeV independent of tan/3, which 
corresponds to a reach in terms of gluino mass of m g ~ 575(750) GeV. 
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Figure 15. Left : The reach of Fermilab Tevatron in the mo vs. mi/2 parameter plane 
of the mSUGRA model, with tan/3 = 10, Ao = and fi > 0, assuming a 5cr signal at 
10 fb — 1 (solid) and a 3a signal with 25 fb^ 1 of integrated luminosity (dashed). The 
red (magenta) region is excluded by theoretical (experimental) constraints. The region 
below the magenta contour has m^ < 114.1 GeV, in violation of Higgs mass limits from 
LEP2. Right: The reach of Fermilab Tevatron in the mo vs. mi/2 parameter plane of 
the mSUGRA model, with tan/3 = 52, Ao = and fi > 0. The red (magenta) region 
is excluded by theoretical (experimental) constraints. The region below the magenta 
contour has < 114.1 GeV, in violation of Higgs mass limits from LEP2. 



4.1.6. Astrophysical and Collider Dark Matter 

Above we have analyzed constraints placed on supersymmetric particle 
physics models, in particular MSSM, by WMAP/CMB astrophysical data. 
The analysis made the assumption that neutralinos constitute exclusively 
the astrophysical DM. It would be desirable to inverse the logic and ask the 
question 42 : "are neutralinos produced at the LHC the particles making up 
the astronomically observed dark matter?" 

To answer this question, let us hrst recall the relevant neutralino in- 
teractions (within the mSUGRA framework) that could take place in the 
Early universe (fig. 16). As we have discussed previously, the WMAP3 




Figure 16. The Feynman diagrams for annihilation of neutralino dark matter in the 
early universe. The Boltzmann factor e~ & M / 20 i n the stau-co-annihilation graph is 
explicitly indicated. 
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constraint (3.8) limits the parameter space to three main regions arising 
from the above diagrams (there is also a small "bulk" region): (1) The 
stau-neutralino (fi — \i) co-annihilation region. Here too is small and 
m i/2 < 1-5 TeV. (2) The focus region where the neutralino has a large 
Higgsino component. Here TO1/2 is small and too > 1 TeV. (3) The funnel 
region where annihilation proceeds through heavy Higgs bosons which have 
become relatively light. Here both too and TO!/ 2 are large. A key element 
in the co-annihilation region is the Boltzmann factor from the annihilation 
in the early universe at kT ~ 20 GeV: exp[-AM/20], AM = M fl - M^o 
implying that significant co-annihilation occurs provided AM < 20 GcV. 

The accelerator constraints further restrict the parameter space and if 
the muon g^-2 anomaly maintains , (c.f. (4.13)), then ji > is preferred 
and there remains mainly the co-annihilation region (c.f. figure 17). Note 
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Figure 17. Allowed parameter space in mSUGRA. Dashed vertical lines are possible 
Higgs masses (from [42]). 



the cosmologically allowed narrow co-annihilation band, due to the Boltz- 
mann factor for AM = 5—15 GeV, corresponding to the allowed WMAP 
range for Vl^ati 2 . 

One may ask, then, whether: (i) such a small stau-neutralino mass 
difference (5-15 GeV) arise in mSUGRA, since one would naturally ex- 
pect these SUSY particles to be hundreds of GeV apart and (ii) such a 
small mass difference be measured at the LHC. If the answer to both these 
questions is in the affirmative, then the observation of such a small mass 
difference would be a strong indication 42 that the neutralino is the astro- 
nomical DM particle, since it is the cosmological constraint on the amount 
of DM that forces the near mass degeneracy with the stau, and it is the 
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accelerator constraints that suggest that the co-annihilation region is the 
allowed region. 

As far as question (i) is concerned, one observes the following: In the 
mSUGRA models, at GUT scale we expect no degeneracies, the AM is 
large, since rriiw governs the gaugino masses, while Too the slepton masses. 
However, at the electroweak scale (EWS), the Rcnormalization Group 
Equation can modify this: e.g. the lightest selectron e c at EWS has mass 
m| c = Too + 0.15m^ 2 + (37GcV) 2 while the \i has mass m| = 
0.16m^ 2 The numerical accident that coefficients of m 2 ^ 2 is nearly the 
same for both cases allows a near degeneracy: for too = 0, e c and \\ be- 
come degenerate at m 1 / 2 — (370-400) GeV. For larger m 1 / 2 , near degeneracy 
is maintained by increasing to to get the narrow corridor in TO0-TO1/2 plane. 
Actually the case of the stau f i is more complicated 42 : large t-quark mass 
causes left-right mixing in the stau mass matrix and this results in the f\ 
being the lightest slepton and not the selectron. However, a result similar 
to the above occurs, with a f\ — Xi co-annihilation corridor appearing. 

We note that the above results depend only on the U(l) gauge group 
and so co-annihilation can occur even if there were non-universal scalar 
mass soft-breaking or non- universal gaugino mass soft breaking at Mq. 
Thus, co-annihilation can occur in a wide class of SUGRA models, not 
just in mSUGRA. Hence, in such models one has naturally near degenerate 
neutralino-staus, and hence the answer to question (i) above is affirmative. 

Now we come to the second important question (ii), namely, whether 
LHC measurements have the capability of asserting that the neutralino (if 
discovered) is the astrophysical DM. To this end we note that, in LHC, the 




Figure 18. SUSY production of neutralinos and decay channels 

major SUSY production processes of neutralinos are interactions of gluinos 
(<?) and squarks (q) (c.f. figure 18), e.g., p + p — > g + q. These then decay 
into lighter SUSY particles. The final states involve two neutralinos x? 
giving rise to missing transverse energy E^ iss ) and four r's, two from the 
g and two from the q decay chain for the example of fig. 18. 

In the co-annihilation region, two of the taus have a high energy ("hard" 
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taus) coming from the \2 ~ ¥ decay (since M^o ~ 2Mf 1 ), while the other 
two are low energy particles ("soft" taus) coming from the f% — * t + Xi 
decay, since AM is small. 

The signal is thus £"™ ISS + jets +r's, which should be observable at the 
LHC detectors 42 . As seen above, we expect two pairs of taus, each pair 
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Figure 19. Number of tau pairs as a function of invariant tt mass. The difference 
Nqs-Nls cancels for mass > 100 GeV eliminating background events (from [42]). 



containing one soft and one hard tau from each x° decay. Since \2 ^ s 
neutral, each pair should be of opposite sign. This distinguishes them from 
SM- and SUSY-backgrounds jets-faking taus, which will have equal number 
of like-sign as opposite-sign events 42 . Thus, one can suppress backgrounds 
statistically by considering the number of opposite-sign events Nos minus 
the like-sign events Nl$ (figure 19). 

The four r final state has the smallest background but the acceptance 
and efficiency for reconstructing all four taus is low. Thus to implement 
the above ideas we consider here the three r final state of which two are 
hard and one is soft. There are two important features: First, Nos-ls in- 
creases with AM (since the r acceptance increases) and Nos-ls decreases 
with Mg (since the production cross section of gluinos and squarks decrease 
with Mg). Second, one sees that Nos-ls forms a peaked distribution. The 
di-tau peak position M^ ak increases with both AM and Mg. This allows 
us to use the two observables Nos-ls and MP° ak to determine both AM 
and M g (c.f. figure 20). As becomes evident from the analysis 42 (c.f. fig. 
21) it is possible to simultaneously determine AM and the gluino mass 
M g . Moreover, one sees that at LHC even with 10 fb _1 (which should be 
available at the LHC after about two years running) one could determine 
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Figure 20. Nqs—LS as function of AM (left graph) and as a function of M g (right 
graph). The central black line assumes a 1% fake rate, the shaded area representing the 
20% error in the fake rate (from [42]). 

AM to within 22%, which should be sufficient to know whether one is 
in the SUGRA co-annihilation region. The above analysis was within the 




Figure 21. Left : Simultaneous determination of AM and Mg. The three lines plot 
constant Nqs_ls and M£° ak (central value and lcr deviation) in the Mg-AM plane for 
the benchmark point of AM =9 GeV and Af§=850 GeV assuming 30 fb luminosity. 
Right: Uncertainty in the determination of AM and Mg as a function of luminosity 
(from [42]). 



mSUGRA model, however similar analyses for other SUGRA models can be 
made, provided the production of neutralinos is not suppressed. In fact, the 
determination of Mg depends on mSUGRA universality of gaugino masses 
at GUT scale, Mq, to relate M^o to Mg thus a model independent method 
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of determining M g would allow one to to test the question of gaugino uni- 
versality. However, it may not be easy to directly measure Mg at the LHC 
for high tan j3 in the co-annihilation region due to the large number of low 
energy taus, and the ILC would require a very high energy option to see 
the gluino. 

One can also measure 42 AM using the signal E^ lss + 2 jets+2r. This 
signal has higher acceptance but larger backgrounds. With 10 fb _1 one can 
measure AM with 18% error at the benchmark point assuming a separate 
measurement of Mg with 5% error has been made. While the benchmark 
point has been fixed in 42 at M g = 850 GeV(i.e. mi/ 2 =360 GeV), higher 
gluino mass would require more luminosity to see the signal. One finds 
that with 100 fb _1 one can probe m 1 / 2 at the LHC up to ~ 700 GeV 
(i.e., Mg up to ~ 1.6 TeV). Finally it should be mentioned that measure- 
ments of AM at the ILC could be made if a very forward calorimeter is 
implemented to reduce the two 7 background. In such a case, AM can 
be determined with 10% error at the benchmark point, thereby implying 
that 42 in the co-annihilation region, the determination of AM at the LHC 
is not significantly worse than at the ILC. 
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Figure 22. Accuracy of WMAP (horizontal green shaded region), LHC (outer red rect- 
angle) and ILC (inner blue rectangle) in determining M x , the mass of the lightest 
neutralino, and its relic density Q x h 2 . The yellow dot denotes the actual values of 
M x and Vt x h? for a sample point in parameter space of mSUGRA: mo = 57 GeV, 
m i/2 = 250 GeV, Aq = 0, tan/3 = 10 and sign(/x) = +1 (from A. Birkedal et ai, 
arXiv:hep-ph/0507214) 

The results on the accuracy of determining DM mass in astrophysics 
and colliders within the mSUGRA framework is given in figure 22. We see 
that the cosmological measurement are at present the most accurate one, 
however, the reader should bear in mind the model-dependence of all these 
results. We now come to demonstrate this point by repeating the analysis 
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for some class of stringy models. 

4.2. Stringy Models and Particle Physics Constraints 

String theory (at least as we know it at present) lives in higher than 
four space-time dimensions (supcrsymmetric strings D=10). Low-energy 
field theory includes gravitational string mutliplet fields (graviton, dilaton 
S(scalar), and supersymmetric partners in Supergravity theories). There 
is an obvious need for compactification to four dimensions. This happens 
dynamically through the Moduli scalar fields, T l , of the string multiplcts, 
which depend on the extra dimensions. 

Originally, it was thought that the requirement of the absence of insta- 
bilities in the vacuum necessitates target-space supersymmetric strings (su- 
perstrings) in general. However, tachyonic instabilities may be welcome in 
cosmological scenarios, hence non supersymmetric target-space-time strings 
may be at play. Target-space supersymmetry needs breaking and must be 
phcnomcnologically consistent, i.e. partners must have masses above a few 
TeV. Consistent breaking of SUSY in string-inspired SUGRA via gaugino 
is possible, and rigorous, and phenomenologically realistic models, in this 
respect, do exist. It is the purpose of this subsection to discuss particle 
physics constraints in the framework of one class of such models. There are 
Modified Constraints on such string-inspired SUSY models from accelerator 
physics and Dark Matter which we shall discuss below. 

The model we shall concentrate upon is a Heterotic string with Orb- 
ifold compactification from ten to four-dimensions with standard model 
gauge group 51/(3) x SU (2) x U(l), three generations, and consistent SUSY 
breaking via gaugino condensate 9 . Below we shall briefly review its most 
important features. 

There is dominance of one-loop soft-SUSY-breaking non-universal 
terms, as a result of superconformal anomalies (non-zero /?- function) . This 
modifies predictions from mSUGRA, and implies and interpolation be- 
tween Anomaly-Mediated SUSY Breaking models and mSUGRA. One dis- 
tinguishes two regimes/scenaria for the SUSY-breaking terms: 
Moduli- dominated: SUSY breaking is driven by the compactification moduli 
fields T l , whose vacuum expectation value (v.e.v.) < T l determines 
the size of the compact manifold . In this regime, there are light scalars 
and relatively heavy gauginos, whose nature depends completely on the 
value of the group-independent coefficient of the universal Green-Schwarz 
counterterm, <5gs 9 - 
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Dilaton- dominated: The dilaton S acquires av.e.v. < S >^ 0, which in 
turn determines the value of string coupling g s at the string scale, and 
transmits, via the (SUSY) auxiliary fields, SUSY Breaking. There are non- 
perturbative corrections to Kahler potential which stabilize the dilaton in 
the presence of gaugino condensation . The associated phenomenology is 
completely different from the moduli-dominated case. We are in a domain 
of heavy squarks and sleptons (of order of the gravitino scale) and light 
gaugino masses, driven by the dilaton auxiliary field v.e.v.'s 




Figure 23. Contours of relative running gaugino masses M3/M2 in the ((Re t) , 5qs) 
plane. These soft masses are at the initial (GUT) scale. The heavy (dark) contour 
is the limit of vanishing gluino mass (there is another such contour in the upper left 
corner on the other side of the self-dual point) . For (Re t) > 1 we also give contours of 
|M 3 /M 2 | = 0.33 (dashed) and 0.75 (solid). 



There is a diverse origin of SUSY Breaking terms in this class of models: 
(i) Some come from the superconformal anomalies , and hence are non- 
universal (proportional to the (3- function of the SU{3) x SU(2) x U(l) 
groups); (ii) some are independent of the gauge group considered (Green- 
Schwarz countcrtcrm, v.c.v. of the condensate). This interplay between 
universality and non-universality implies a rich phenomenology, new trends 
in the search of supersymmctric particles in accelerator and astro-particle 
physics. 

The heterotic string models we shall analyze here make use of the 
so-called Pauli-Villars -Anomaly Mediated Supersymmetry Breaking (PV- 
AMSB) scenario 9 . There is an important feature of the AMSB-string mod- 
els which turns out to be important for DM searches, namely the presence 
of non-thermal LSP relic densities. Indeed, in AMSB studies it was found 
that the neutralino thermal relic density is generically too small to explain 
the amount of dark matter: the wino content of the LSP is quite high. 
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tanp=5 tanp = 35 




b + b + 

Figure 24. Upper: Constraints on the moduli-dominated parameter space for tan/3 = 5 
(left) and tan /3 = 35 (right) with p = and (Re t) =2.0. Constraints on the (M3/2, &gs) 
plane are given for fi > 0. The dark shaded regions on the left have a stau LSP. The 
tan /3 = 35 plot also has a region with a gluino LSP. For tan /3 = 5 the region labeled 
"O" has the cosmologically preferred relic density of neutralinos. No such region exists 
for the higher tan /3 plot. In that case the exclusion contours are due to (from bottom 
right to upper left) Charge Conjugation Breaking vacua, the chargino mass, too large 
SUSY contributions to (g M — 2), the Higgs mass limit and too large a b — > S7 rate. 
Lowcr: Constraints on the dilaton-dominated parameter space for tan /3 = 5 (left) and 
tan/3 = 35 (right). Constraints on the &+) plane are given for fj, > [with 6+ 

the largest /3-function coefficient among condensing gauge groups of the hidden sector] . 

Additionally, co-annihilation between the LSP and the lightest chargino is 
also very efficient. Both of these effects combine to make the thermal relic 
density of LSP negligible. Thus the anomaly-mediated character of the 
gaugino sector in this model necessitates a non-thermal production mech- 
anism for neutralino LSPs, or another candidate for the cold dark matter 
must be postulated. 

We next remark that, in the moduli-domination scenario of the heterotic 
string models, the various Kahlcr moduli are not stabilized, i.e. (Re t) is not 
fixed, and one must use the soft terms. The value of the Green- Schwarz 
coefficient 5qs, then, becomes relevant for the determination of gaugino 
masses 9 , which are "running" with it (c.f. figure 23). 

The relevant astro-particle physics constraints for one indicative exam- 



pie of hctcrotic orbifold models 9 , in some range of the parameter space, 
are summarized in figure 24. The reader's attention is drawn to the strong 
suppression (or disappearance) of thermal neutralino relic densities in some 
regions of the parameter space. 

4.3. Non-critical (non-equilibrium) Stringy Q- Cosmology 

The results described on the previous sub section were based on a critical 
string theory, where the dilaton field is stabilized. However, one may en- 
counter situations in some non-equilibrium stringy cosmologies, described 
by non-critical strings 10 (Q-cosmologies), in which the dilaton is not sta- 
bilized. Such cosmologies might arise, for instance, in colliding brane world 
scenarios, where the cosmically catastrophic early universe event of colli- 
sion suffices to induce departure of the associated string theory, describing 
various excitations on our brane world, from conformal invariance (on a 
world sheet) and hence non-criticality. In such models, the dilatons are 
time dependent and are not stabilized. At late eras of the Universe, such 
a time dependent dilaton may behave as quintessence-like field leading to 
acceleration of the Universe 43 ' 10 > 17 . The departure of criticality has also 
other consequences for the low-energy field theory, namely the existence of 
off-shell terms 10 , i.e. the variations of the relevant effective action with 
respect to the fields in the gravitational string multiplet (gravitons, dila- 
tons) may be non zero. A detailed analysis 44 shows that in such models, 
where the dilaton is not stabilized, it is possible to have thermal relic abun- 
dances of DM particles, such as neutralinos in supersymmetric cases, whose 
density, however, obeys a modified Boltzmann equation by the dilaton and 
off-shell source terms: 

f + 3 (£) n = $n - \ (e- + 2e*/3*) n + J ^C[f] (4.14) 

with n = J d 3 p / the species number density. The /3-terms denote the off- 
shell terms in the gravitational (graviton, dilaton) multiplet of the string: 
/vJ av « ap^j h fif 7^ 0, where S is the low-energy string-inspired 
target-space effective action. 

If we set, for brevity, T(t) = $ - \ (e~* '^ rav + 2e*/3 <I, ) , acting as a 
source term on the right-hand-side of the Boltzmann equation (4.14), then 
we can solve this equation to obtain for the thermal DM relic abundance 
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(assuming a single species, e.g. neutralino \) 44 : 

n*h 2 = (^hl) no _ source x (|) V2 exp Qf ' (4-15) 

with (fi^/ip) no _ source the standard (equilibrium) cosmology result (3.10); 
the star notation denotes quantities at the freeze-out point, g = g + 
!§-T -4 Ap, with T the temperature, g is the effective number of (ther- 
mal) degrees of freedom in this non-critical string Universe and Ap de- 
notes collectively the dilaton and off-shell (due to the /3-functions) terms 
in the effective modified Friedman equation of the Q-cosmology 10 > 44 : 
H 2 = (p + Ap). 

Depending on their signature, the Source Terms T have different influ- 
ence on the relic abundance, with profound consequences on the prospects 
for detecting supersymmetry in such models at colliders. For instance, for 
a given model 44 and within a certain region of the parameters, there is 
a reduction of neutralino relic abundance, as compared with conventional 
cosmologies, by a factor of about 1/10 (c.f. figure 25), leading to a relax- 
ation on some of the constraints regarding SUSY detection prospects at 
LHC imposed by mSUGRA models. 



■uj. 
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Figure 25. Left: In the thin green (grey) stripe the neutralino relic density is within 
the WMAP3 limits for values Ao = and tanfi = 10, according to the source-free T=0 
conventional Cosmology. The dashed lines (in red) are the lcr boundaries for the allowed 
region by the g — 2 muon's data as shown in the figure. The dotted lines (in red) delineate 
the same boundaries at the 2fr's level. In the hatched region 0.0950 > n C DMh 2 , while 
in the dark (red) region at the bottom the LSP is a stau. Right: The same as in left 
panel, but according to the non-critical-string calculation, in which the relic density is 
reduced in the presence of dilaton sources r = tj> ^ 0. 
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5. Conclusions and Looking Ahead 

In this set of lectures I have reviewed various astrophysical methods for con- 
straining the Universe energy budget, and used such results to constrain in- 
teresting particle physics models of cosmological relevance. I have discussed 
the issue of calculating thermal-relic DM abundances in those models, with 
the intention of placing constraints on interesting particle physics models, 
such as supersymmetry, exploiting astrophysical data. I mainly concen- 
trated on DM searches at LHC. I came to the conclusion that LHC could 
shed some light on the issue as to whether the astronomically observed DM 
is the neutralino, but this is a highly model dependent statement. I have 
analyzed briefly the phenomenology of various SUGRA models (mSUGRA, 
and some string inspired ones, including some relaxation non-equilibrium 
dilaton-quintessence models). The associated phenomenologies are very dif- 
ferent, depending crucially on the details of the underlying dynamics, such 
as the type of the SUGRA model, the way of SUSY breaking etc. There are 
model independent methods for testing DM at colliders, but they pertain 
to subdominant processes at colliders. 

For future directions it would be desirable to explore in more detail 
SUSY models with CP violation, which recently started attracting atten- 
tion 45 , since, due to bounds on Higgs huh > 114 GeV, we now know that 
the amount of CP Violation in the Standard Model is not sufficient to gen- 
erate the observed baryon asymmetry of the Universe 46 , and hence SUSY 
CP violation might play an important role in this respect. At this point I 
mention that parameters in SUGRA models that can have CP phases are 
the gaugino and higgsino masses and trilinear sfermion-Higgs couplings. 
CP phases affect co-annihilation scenaria, and hence the associated par- 
ticle physics dark matter searches at colliders 45 . Another direction is to 
constrain SUSY GUTs models (e.g. flipped SU(5)) using astrophysical 
data 7 , after taking, however, proper account of the observed dark energy 
in the Universe. Personally, I believe that this dark energy is due to some 
quintessence (relaxing to zero (non-equilibrium) field). WMAP data point 
towards an equation of state of quintessence type, w = p/p — ► —1 (close to 
that of a cosmological constant, but not quite —1). Such features may be 
shared by dilaton quintessence, as discussed briefly above in the context of 
string theory. The issue is, however, still wide open and constitutes one of 
the pressing future directions for theoretical research in this field. 

On the experimental side, LHC and future (linear) collider, but also 
direct 23 , dark matter searches could shed light on the outstanding issue 
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of the nature of the Cosmological Dark Sector (especially Dark Matter), 
but one has to bear in mind that such searches are highly theoretical- 
model dependent. To such ideas one should also add the models invoking 
Lorentz violation as alternative to dark matter. Clearly, particle physics 
can play an important role in constraining such alternative models in the 
future, especially in view of the upcoming high-precision terrestrial and 
extraterrestrial experiments, such as Auger, Planck mission, high-energy 
neutrino astrophysics etc. 

Nothing is certain, of course, and careful interpretations of possible 
results are essential. Nevertheless, the future looks promising, and certainly 
particle physics and astrophysics will proceed together and provide a fruitful 
and complementary experience to each other and exchange interesting sets 
of ideas for the years to come. 
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